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PREFACE 


This  report  is  published  to  provide  guidance  to  coastal  engineers  in 
planning  wave  data  collection  in  coastal  waters  for  climatology  purposes, 
including  wave  direction.  The  popularity  of  three-gage  arrays  for  pro- 
posed wave  direction  measuring  systems  makes  it  necessary  to  evaluate 
the  capabilities  and  the  limitations  of  these  arrays.  The  availability 
of  the  CERC  five-gage  array  at  Pt.  Mugu,  California,  provided  a unique 
opportunity  for  evaluating  the  performance  of  wave  recording  systems  and 
the  directional  capabilities  of  three-gage  arrays.  The  work  was  carried 
out  under  the  wave  measurement  and  analyses  program  of  the  U.S.  Army 
Coastal  Engineering  Research  Center  (CERC) . 

Tliis  report  was  prepared  by  Dr.  Dinorah  C.  Esteva  under  the  super- 
vision of  Dr.  D.  Lee  Harris,  Chief,  Coastal  Oceanography  Branch,  Research 
Division.  The  author  acknowledges  the  valuable  insight  and  comments  pro- 
vided by  Dr.  D.  Lee  Harris  and  Mr,  R.  P.  Savage,  Chief,  Research  Division, 
CERC. 

Comments  on  this  publication  are  invited. 
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COm'ERSION 

U.S.  customary  units 
to  metric  (SI)  units 

FACTORS,  U.S. 
UNITS  OF 

of  measurement 
as  follows: 

CUSTOMARY  TO  METRIC  (SI) 

MEASUREMENT 

used  in  this  report  can  be  converted 

Multiply 

by 

To  obtain 

inches 

25.4 

mi llimeters 

2.54 

centimeters 

square  inches 

6.452 

square  centimeters 

cubic  inches 

16.39 

cubic  centimeters 

feet 

30.39 

centimeters 

0.3048 

meters 

square  feet 

0.0929 

square  meters 

cubic  feet 

0.0283 

cubic  meters 

yards 

0.9144 

meters 

square  yards 

0.836 

square  meters 

cubic  yards 

0.7646 

cubic  meters 

miles 

1.6093 

kilometers 

square  miles 

259.0 

hectares 

knots 

1.8532 

kilometers  per  hour 

acres 

0.4047 

hectares 

foot-pounds 

1.3558 

newton  meters 

millibars 

1.0197  X 

10 

kilograms  per  square  centimeter 

ounces 

28.35 

grams 

pounds 

453.6 

grams 

0.4536 

kilograms 

ton,  long 

1.0160 

metric  tons 

ton,  short 

0.9072 

metric  tons 

degrees  (angle) 

0.1745 

radians 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  Kelvins^ 

^To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  readings, 
use  formula:  C = (5/9)  (F  -32). 

To  obtain  Kelvin  (K)  readings,  use  formula:  K = (S/9)  (F  -32)  + 273.15. 
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HVALUATION  OF  THE  COMPUTATION  OF  WAVE  DIRECTION  WITH  THREE-CAGE  ARRAYS 

by 

Dinorah  C.  Esteva 
I.  INTRODUCTION 

Wave  direction  is  an  important  parameter  in  the  solution  of  many 
coastal  engineering  problems.  A knowledge  of  wave  direction  is  essen- 
tial for  (a)  estimating  the  direction  and  magnitude  of  sediment  trans- 
port by  waves,  (b)  using  refraction  calculations  to  infer  wave  conditions 
at  one  site  from  measurements  made  elsewhere,  and  (c)  verifying  theories 
of  wave  generation. 

Visual  observations  of  wave  directions  have  been  collected  by  ship- 
board observers  for  over  a century.  About  20  years  ago  the  Beach  Ero- 
sion Board  (BEB) , predecessor  to  the  Coastal  Engineering  Research  Center 
(CERC) , engaged  the  assistance  of  U.S.  Coast  Guard  installations  in  the 
collection  of  visual  observations  of  breaker  direction  from  shore.  How- 
ever, objective  determinations  of  wave  direction  are  desirable  without 
being  restricted  to  location,  time  of  day,  or  visibility  condition.  The 
capability  to  do  so  involves  the  use  of  wave  measuring  instruments. 

Panicker  (1971,  1974)  presents  extensive  reviews  of  reports  dealing  with 
the  determination  of  wave  direction  from  instrument  records  with  partic- 
ular emphasis  on  those  involving  sea-surface  elevation  or  pressure  records. 

In  March  1970,  CERC  installed  an  array  of  wave  gages  at  Pt.  Mugu, 
California.  Records  from  the  array  were  to  be  used  to  compare  redundant 
values  of  wave  direction  and  to  estimate  the  level  of  accuracy  of  the 
computations.  The  available  procedures  for  determining  wave  direction 
from  an  array  involved  assumptions  that  had  not  been  thoroughly  estab- 
lished. Tlius,  the  records  from  the  array  would  also  be  used  in  a syste- 
matic examination  of  these  assumptions,  and  of  the  reliability  of  wave 
gages . 

This  study  discusses  the  array  performance  and  the  information  gained 
about  wave  direction.  Redundant  values  of  directions  were  obtained  from 
the  10  three-gage  arrays  possible  with  five  gages.  The  mathematical  model 
used  assumes  that  the  sea  surface  is  the  result  of  the  superposition  of 
a small  number  of  narrow-banded  wave  trains  consisting  of  long-crested 
waves  traveling  in  well-defined  directions.  It  was  also  assumed  that 
only  one  wave  train  is  present  with  a particular  period.  The  first 
assumption  is  supported  by  the  energy  spectra  computed  at  CERC  (Thompson, 

1974)  , by  aerial  photos  of  the  sea  surface  (Fig.  1) , and  by  radar  images 
of  the  wave  field  (Fig.  2).  Many  published  reports  include  photos  similar 
to  that  in  Figure  1;  e.g.,  McClenan  and  Harris  (1975).  Fujinawa  (1974, 

1975)  conjectured  that  narrow  directional  spread  might  be  responsible 
for  the  incomplete  recovery  of  the  true  directional  spectrum  from  field 
records  in  his  computations  using  high  directional  resolution. 

Average  values  of  wave  direction  for  bands  of  constant  frequency 
width  were  computed  from  cross-spectra  between  gage  pairs.  Direction 
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Figure  2.  Radar  scan  of  wav'C  field  at  Xausct,  Massachusetts. 


estimates  for  all  hands  0.01  hert  a wide  between  approximate])’  30  and  3 
seconds  were  obtained  for  the  10  arra>'s.  The  results  displ.i)’cd  discrep- 
ancies of  the  order  of  d0°  for  those  hands  with  central  periods  ahov'e  KJ 
seconds  and  of  180°  for  those  with  shorter  central  periods.  It  had  been 
expected  that  the  array  would  yield  direction  to  better  than  20°  and  for 
periods  between  25  and  7 seconds. 

To  isolate  problems  associated  with  the  cal 'ulations , tlic  propagation 
of  narrow-banded  wave  trains  across  the  array  was  simulated  in  a computer. 
ITie  computational  model  was  applied  to  the  simulated  observations  using 
tlie  maximum  fi'  quency  resolution  available  from  spectral  computations 
based  on  20-minute  records.  It  was  found  that  the  directional  results 
obtained  with  this  model  are  highly  dependent  on  the  spectral  width, 
both  in  frequency  and  direction,  of  the  wave  train  involved  and  on  the 
relationship  between  wavelength  at  the  site  and  gage  separations.  Tlie 
assigned  directions  were  recovered  witliin  1°  for  16-second  waves  when  the 
frequencies  of  the  spectral  components  in  the  wave  train  differed  by  O.'^O.i 
licrtz  or  more,  and  the  directions  were  spread  witliin  a 5°  arc.  This  fre- 
quency separation  results  in  a minimum  difference  in  periods  of  0.7  second 
for  waves  with  periods  near  16  seconds. 

Application  of  the  same  analyses  to  field  wave  pressure  records  with, 
standard  deviations  above  0.61  meter  (2  feet)  resulted  in  an  average  dis- 
crepancy of  20°  among  computed  directions  for  narrow-banded  wave  trains 
with  periods  longer  than  10  seconds.  Larger  discrepancies  resulted  for 
shorter  periods.  Thus,  accuracies  no  better  than  20°  can  be  expected  for 
wave  directions  resulting  from  three-gage  arrays. 

1 . The  System. 

A minimum  of  three  gages  is  required  for  a unique  determination  of 
wave  direction  by  most  proposed  models.  Since  these  models  make  a few 
assumptions  about  the  nature  of  ocean  waves  which  have  not  been  estab- 
lished, some  redundance  was  thought  to  be  necessary  whicli  would  require 
a minimum  of  four  gages.  However,  it  was  agreed  that  a five-gage  array 
would  increase  the  probability  of  redundance  in  the  ocean  environment. 

An  array  wa;.  designed  at  CLRC’  by  Leon  E.  Bergman,  s t at  is  t i c L an- engineer  , 
while  o::  sabiiatical  leave  from  the  Universit)-  of  California,  Borkelo)’, 
when  the  experiment  was  being  planned,  lie  investigated  the  di  rcct  iciia.  1 
resolving  jidwer  of  several  arr.’iy  geometries  and  concluded  tiiat  the  jiat- 
tem  shown  in  I'igure  .3  would  lie  tlie  most  suitable  for  the  conditions  to 
lie  exjiccted  at  i’t.  Mugu  (Borgm.’ui  .ind  I’anicker,  1P70). 

The  arr.’i)'  was  installed  off  I’t.  Mugu,  approximately  80.47  kilometers 
[50  miles)  northwest  of  Los  Angeles  (Tig.  4),  in  about  0.14  meters  i 50 
feet)  of  water,  0.7(i  meter  (2.5  feet)  from  the  bottom.  Tlie  gages  in  the 
arm)'  are  pressure  t r.ansduccrs  developed  mostly  at  CliRC  (Williams,  lOh'.!}. 
The  heart  of  the  system  is  a Fairchild  pressure  transducer  ii-hicii  is  potted 
inside  a 2-inch  Plexiglas  tulu-  (Fig.  5)  (Peacock,  1074).  A plastic  tube 
filled  with  silicone  oil  transmits  the  pressure  fi'om  the  seawater  to  tlie 
pressure  transducer.  llu'  silicone  oil  is  separated  from  seawater  by 
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rubber  diaphragms.  A Teflon  scouring  pad  dipped  in  ;inti  fouling  paint 
separates  the  rubber  diaphragm  from  the  end  of  the  package  which  admits 
the  seawater  pressure.  The  instrument  is  mounted  vertically  on  a tri- 
pod (Fig.  6)  . Tl\e  signals  from  the  pressure  transducers  are  brought  by 
cables  to  a recording  and  transmitting  station  onshore. 

2.  Data  Collection. 


The  array  went  into  operation  on  27  Marcli  1970.  Tl>c  water  pressure 
at  the  five  gages  was  registered  continuously  at  a rate  of  four  times  a 
second  on  digital  magnetic  tape.  Difficulties  experienced  during  the 
first  year  with  individual  sensors  were  mostly  of  short  duration  and 
were  presumed  to  be  due  to  biological  activity.  However,  major  diffi- 
culties were  experienced  with  the  recording  system,  and  on  16  Marcli  1971 
the  recorder  was  disconnected  at  Pt.  Mugu  and  transferred  to  the  CFRC 
laboratory,  then  located  in  Washington,  D.C.  Recording  from  all  five 
gages  was  reinitiated  at  CERC  on  9 April  1971.  Records  were  obtained 
continuously  until  3 January  1972  when  the  recorder  failed.  During  this 
period  much  of  the  data  were  useless  because  of  an  unacceptable  level  of 
high-frequency  noise.  Ttie  source  of  the  noise  was  difficult  to  locate 
and  was  not  eliminated  from  the  signal  until  shortly  before  the  failure 
of  the  recorder  in  January.  Since  2 February  1972,  records  from  throe 
of  the  five  gages  at  Pt.  Mugu  have  been  included  in  the  time-shared 
recording  of  waves  from  east  and  gulf  coast  wave  station;;  (Peacock, 
1974).  In  1972,  data  were  recorded  for  20  minutes  out  of  each  hour; 
since  February  1973,  data  have  been  recorded  for  20  minutes  out  if  each 
2-hour  interval. 


II.  FIELD  DATA  .ANALYSIS 

The  five  gages  in  the  array  provided  uninterrupted  data  for  most  of 
the  first  year  of  operation.  Eight  daily  observations,  each  consisting 
of  simultaneous  20-minute  records  from  the  five  gages,  ivierc  processed 
from  these  data.  The  observations  had  starting  times  wit'un  1.5  hours 
of  the  weather  synoptic  times  (0100,  0400,  0700,  1000,  1300,  1600,  and 
1900  hours  P.s.t.).  The  potential  energy  in  the  wave  field  is  propor- 
tional to  the  variance  of  the  time  history  of  sea-surface  elevation  at 
a fixed  location  (Kinsman,  1965).  For  most  conditions,  the  standard 
deviation  of  the  surface  displacement  is  one- fourth  of  the  significant 
wave  height.  The  standard  deviation  of  the  pressure  at  a fixed  depth  is 
roughly  proportional  to  the  wave  height  and  may  also  serve  as  a measure 
of  the  wave  height. 

The  standard  deviation  of  the  recorded  pressure  was  computed  from  the 
records  of  each  of  the  five  gages  for  eight  records  each  day.  The  stand- 
ard deviation  of  the  record  from  each  gage  was  compared  witli  the  average 
of  the  five  gages.  If  the  standard  deviation  from  any  gage  differed  from 
the  average  by  more  than  20  percent,  the  record  from  that  gage  was  deleted 
and  the  average  recomputed  from  the  remaining  gages.  .A  comparison  of  the 
individual  stand.ird  deviations  from  the  mean  for  that  time  period  is  shown 
in  Table  1. 
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Table  1.  Percent  of  observations  where  the  largest 
departure  of  the  standard  deviations  from 
the  moan  in  the  observations  was  as  indi- 
cated (871  observations  in  1970). 


Deviation  from  mean  (pet) 

Percent  of  Observations 

< 2 

41 

3 to  10 

52 

11  to  20 

1 

> 20 

6 

This  com.parison  indicates  the  system  operated  consistently.  Tlie 
field  wave  records  used  for  wave  direction  computation  (discussed  later) 
were  chosen  from  the  observations  in  Table  1,  for  which  the  standard 
deviations  from  all  gages  differed  by  3 percent  or  less  from  their  mean 
and  for  which  the  average  significant  wave  height  (uncompensated  for 
attenuation  with  depth)  was  above  0.61  meter. 

Fourier  analysis  provides  a reliable  procedure  for  obtaining  the 
periods  of  the  most  important  waves.  A Fourier  analysis  of  sea-surface 
elevation  (or  pressure)  with  time  results  in  the  distribution  of  energy 
with  frequency,  usually  referred  to  as  the  energy  speetrum.  The  energy 
spectrum  for  the  record  from  gage  5 in  each  of  the  eight  daily  observa- 
tions was  computed  using  the  Fast  Fourier  Transform  algorithm  developed 
by  Cooley  and  Tukey  (1967).  The  first  1,024  seconds  of  the  20-rainute 
record  was  used  in  this  computation.  Gage  5 was  chosen  because  of  a 
good  history  of  performance. 
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Fast  Fourier  Transform  computations  yield  the  contribution  to  the 
variance  at  each  of  a set  of  frequencies  which  are  harmonics  of  a funda- 
nental  given  by  the  inverse  of  the  record  duration,  T.  In  this  study, 
the  frequencies  of  these  harmonics  are  referred  to  as  speatral  freauenaies , 
ind  the  corresponding  periods,  given  by  their  inverse,  as  speatral  periods. 
’he  energy  content  between  32  seconds  and  3 seconds  was  used  to  normalize 
he  spectrum.  The  lower  limit  on  period  was  estimated  from  the  thickness 
>f  the  water  column  above  the  pressure  sensors.  A summarized  spectrum 
Fig.  7)  was  formed  by  combining  the  energy  content  in  11  adjacent  spec- 
ral  periods.  The  band  width  in  the  summarized  spectrum  was  slightly 
arger  than  10”^  hertz  (0.0107  hertz).  The  energy  appearing  at  each 
pectral  period  in  the  pressure  spectrum  was  compensated  for  attenuation 
ith  depth  by  using  the  classical  hydrodynamic  pressure  correction: 


F(k,h) 


cosh  kh 
cosh  kAz  ’ 


(1) 


here  k is  the  wave  number,  h the  water  depth,  and  Az  the  vertical 
istance  of  sensor  from  bottom.  This  resulted  in  a surface  or  compensated 
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.spectrum.  Surface  root-mcan-squarc  (rns)  values  were  obtained  from  tlie 
compensated  energy  at  each  spectral  period. 


Since  1972,  four  simultaneous  20-minutc  records  from  eacli  of  tlie 
tliree  gages  included  in  t)ic  time-shared  recording  sec|uetuc'  are  processed 
daily.  More  records  arc  processed  during  special  studies.  The  records 
processed  from  gages  1,  2,  and  7>  of  tlie  Pt.  Mugu  arra>'  start  al)out  OKiO, 
0700,  1300,  and  1900  hours  (P.s.t.).  The  significant  wave  lieight,  tlie 
distribution  function,  tlie  first  five  moments  of  the  distribution,  and 
the  pressure  spectra  are  computed  for  tliese  records  in  a study  of  wave 
record  variability.  These  records  are  not  analysed  for  wave  direction. 

1 . Computation  of  Wave  Direction. 

For  a long-crested  sinusoidal  wave  with  frequency  jiropagating 

in  direction  (Fig.  8) , the  phase  difference  between  locations  1 and 

2 with  coordinates  (xj,yi^)  and  (X;,,y2)  risspectively , is  given  by; 

”^12  = (>'l  ■ >'2) 

where  = 2r/L^  is  the  wave  number  associated  with  frequency  a^, 

and  is  the  wavelength.  The  subscript  is  used  to  indicate  the 

possible  presence  of  different  wave  trains  with  different  frequencies 
and  directions. 


The  addition  of  the  wave  profile,  03,  at  a third  noncol inear  loca- 
tion allows  solving  for  the  sine  and  cosine  of  Thus,  a unique  solu- 

tion for  the  wave  direction  is  obtained  from  the  following  equation  when 
the  signs  of  numerator  and  denominator  are  considered.  (see  App.  A.); 


tan"  ^ 


[(Xi  - X3)  $33 

[(Tl  - >'2) 


(X3  - X3)  $33]/D 

(y\  - 


(3) 


where  $33  is  the  phase  difference  between  the  third  and  first  locations, 
and  D is  a function  of  gage  separation. 

Phase  differences  between  locations  for  each  different  wave  period 
are  the  only  unknowns  in  the  right-hand  side  of  equation  (3).  Estimates 
of  a representative  phase  difference  between  gage  pairs  for  bands  of 
constant  frequency  width  are  easily  computed  from  cross-spectra  of  the 
wave  (pressure)  records.  These  spectra  give,  for  each  band,  average 
values  of  the  covariance  of  the  wave  records  along  two  perpendicular 
direct i ons . 


Substitution  of  those  representative  values  of  phase  difference  into 
equation  (3)  affords  an  expedient  and  economic  means  of  obtaining  esti- 
mates of  a "representative"  wave  direction  for  each  of  the  spectral 
bands,  provided  the  results  arc  of  engineering  ase.  The  agreement  among 
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Figure  8.  Long- crested  wave  propagating  in  direction 


redundant  computations  of  direction  from  tlie  10  three-gage  arrays  is  an 
indication  of  the  degree  of  confidence  that  can  be  placed  on  the  result- 
ing directions. 


Directions  were  computed  for  a few  observations  (using  the  equation 
below)  for  the  representative  phase  difference  in  each  spectral  band, 
where  again,  the  sign  of  numerator  and  denominator  mast  be  considered. 


12," 


tan 


- 1 


(4) 


The  subscript  A in  equation  (4)  indicates  t!ic  value  is  for  a given 
spectral  band.  The  cospectrum,  Co,  and  quadrature  spectrum,  Quad, 
are  defined  as; 


Co^  = ).  XjX2  cos  $J2 
A 

Quad^^  = f XjX^  sin  * (5) 

A 

where  Xj  and  X2  arc  the  spectral  amplitudes  from  wave  records  1 and 
2,  and  5.  indicates  summation  over  the  adjacent  spectral  periods  com- 
bined to  make  up  a band.  Jenkins  and  Watts  (1908)  showed  that  the  defi- 
nitions in  equation  (5)  are  equivalent  to  the  more  standard  definitions 
based  on  correlation  functions. 

A computer  output  for  runs  using  this  approach  is  presented  in  Appen- 
dix B.  Summaries  of  the  results  of  the  observations  at  0700  hours,  25 
June,  and  2010  hours,  28  June,  are  given  in  Table  2.  The  first  column 
in  the  table  gives  the  period  at  the  center  of  the  band;  the  second 
column  gives  the  average  percent  energy  in  each  band  for  the  five  gages. 
Ihe  last  10  columns  give  the  comfiuted  "representative"  direction  of  each 
band  for  the  three-gage  array.  Directions  are  measured  from  the  seaward 
normal  with  positive  values  counterclockwise  and  negative  values  clock- 
wise. The  table  shows  that  in  these  two  observations  disagreements  in 
direction  of  the  order  of  20°  are  obtained  for  the  longer  period  peaks 
and  of  100°  for  the  shorter  period  peaks.  Hxamination  of  the  computer 
output  in  Appendix  B indicates  the  results  arc  typical. 

The  array  had  been  exjjected  to  yield  directions  to  better  th,in  20° 
for  periods  between  25  and  7 seconds.  Understanding  of  the  problems 
involved  was  souglit  by  simulating  the  propagation  through  the  array  of 
narrow-banded  wave  trains  traveling  in  a specified  direction  (discussed 
in  next  subsection). 

2.  Simulated  Data  Analysis. 

In  simulating  the  wave  records  for  use,  special  consideration  was 
given  to  wave  period  and  to  the  difficulties  arising  in  spectral  analysis. 
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Assume  a sinasoid  with  frequency  given  as: 


d’T  (m  + '') 
NA  t 


(hJ 


where  At  is  the  interval  of  time  l)ctween  samples,  | ; is  less  than 

or  equal  to  1/A,  iind  N'At  = T,  t)>e  record  duration, 

liquation  ((i)  provides  for  assigning  frequencies  sliieh  differ  from  the 
spectral  frequencies.  fhe  contribution  to  the  variance  .it  the  spectral 
frequencies  of  tlie  sampled  record  is  given  by  .S,^  as  ; 


where  and  b.,  are  the  Fourier  coefficients. 


Harris  (1974]  showed  that  for  values  of  m near  m (i.c.,  for  spec- 
tral frequencies  near  tlie  frequency  of  the  sinusoid),  and  for  m far 
removed  from  one  and  N/A,  the  approximations  below  are  good  estimates 
to  the  coefficients. 


a,.,  i 


b„.  = 


A sill  i'  cos  (J.  - ) 

( m - in  + ' ) 

A sin  ■■  sin(i  - : ■ ') 
■•  ( m - m + 1 


(8) 


Slow  convergence  of  the  energy  toward  the  s|-)eetral  ;'eriod  closest  to 
the  assigned  period  is  clearly  indicated.  Thus,  the  energy  is  spread 
over  adjacent  spectral  periods.  Tliis  spreading,  due  to  tiie  finiteness 
of  tlie  record,  is  usually  referred  to  as  e.  A Ac  yi’i'. 

The  technique  routinely  used  at  ChRC  to  decrease  spillo\’cr  is  to 
apply  tlie  cosine  bell  data  window  as  defined  by  : 

1 - cos  y;  , i = 1,  . . . N , (9) 

where  arc  the  values  in  the  original  record.  The  Fourier  coefficient 

for  the  resulting  function  yi,  is  given  by: 

y • .A  s i n Tw  cos(  v - r i] 

‘*'v  “ , , , ’ 

A a ( m - m + \)  [ ( m - m * '■  ] ' - 1 ] 


S/ 

b 


A sin  .1  '■  s i n ( ; - a n ) 

27T(m  - m + 8)  [(n  - m - 1] 


Thus,  convergence  is  greatly  increased  and  spillover  is  effectively 
reduced  to  tliree  adjacent  spectral  periods. 


i 

i 
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Ihe  os  ;no  !s' i 1 data  window  was  ajijilied  to  the  simulated  recoi'ds ; 
tiieia'i'e  .-e , it  was  sufficient  to  consider  wave  trains  consisting  of  three 
sinusoid.-',  with  ne,'irh\'  periods  spread  over  at  most  six  adjacent  s])ectral 
pi't'iods.  In  general,  the  periods  of  waves  in  the  ocean  will  differ  from 
the  spectra]  periods.  Tlnis , tlie  sinusoids  in  each  simulated  ohservation 
were  sjiec  i f i ca  1 1 assigned  periods  differing  slightl)'  from  tlie  spectral 
periods  l'\'  use  of  eipiation  (<))•  Centra]  periods  of  ahout  8 and  1()  seconds 
were  clio:;en  to  sinmlate  8-  and  l(^-second  swells.  Swells  with  periods  in 
tlii,''  range  a -e  oiiserved  on  the  west  coast. 

i .acii  ■■  i i:ul  at  eii  wave  record  consisted  of  values  of  the  wave  profile 
at  t’K-  f i \ e reige  locations  comjiuted  at  0 . .’s- second  intervals  for  17.07 
minute^  1 . M seconds),  to  simulate  the  -.ampling  rale  and  recorii  dura- 
•■ion.  cu  • .man  ly  used  at  CfMC  for  field  d,.ii".. 

; IK  t'i'-.  . siaomids  were  assigue(,l  spc-cific  directions  and  tero  initial 
phase  at  till.'  or,  ;in  af  coordinates,  and  were  propagated  across  the  arra>' 
assui’iin  t .'.mst.Uit  depth  ot  0.14  meters.  .appendix  ti  shows  that  t!ie 
!■  i n : t (.  iourier  Iraiisform  gives  the  I'orrect  phases  for  three  sinusoids 
•i;U.  eo.ifiaed,  jiroi'ided  the  sinusoids  arc  assie.neti  the  same  direction, 

;;r.irl>'  ' iie  saaie  .unplitudes,  aiui  tlte  freqiK-ncy  difference  of  each  coiii- 
pur.eiit  ro  li'  ' nearest,  spoLtiai  frequence  i thi,'  same'.  \ fretpienc)'  di  f- 
ference  -f  ' . 1 51 /' 1 , 0 J 1 f O.ooni.'O  li.-rt.-.  '.as  ■...osc-':. 

' ’la  I'aci  e' r’ s t 1 cs  of  the  first  eigltt  sini'dated  observations  a v.' 
in  la'le  5,  Kmigli  cons  i de  r it  i ons  of  refr.iciie'ii  ;;;ini;  linear  fheor\' 

( llcC  1 i.-nan , )ie'ld  from  tlie  normal  le  tin.  coastline  as  the  inaxi- 

mum  lOossible  direction  that  waves  with  a lo-seconel  pwriod  may  have  at 
tlie  Jeptli  of  the  .iiTa).  Inri'ctiotn  within  71'^  aiul  7 1"  S.  from  tiic* 
normal  iscta.'  chosen  for  the  first  feiur  s i inul  ,it  e-el  wave'  trains.  Since  waves 
witli  an  8-se  ,ionu  perioel  may  approach  the'  ce>as f 1 i ne  at  the  arrax'  'ite'  from 
a much  wide'r  arc,  eiirectietns  up  to  ,N'.  ;iinl  nO"  S.  were  ..ise'd  for  tlie 

directions  of  the-  fi  ftli  to  ei;;htl)  simulated  ti-ain.s. 

The  een.i'puted  sp'.-erra  for  tle-.e  : gb.  t simulated  observat  i enis  are  shown 
in  .'Xpjiendix  1.)  (l-'igs.  D-  1 to  In  these  figures  the  veiriance  of  the 

record,  proportional  to  the  ener;'.}',  at  eaeii  spectra.l  period  is  j'lottol 
versu.s  a lineai'  frequency  scale.  No  greniiting  of  i tu'  variance  at  adjacent 
spectral  periods  lias  Ixeen  made.  These  ungroupe,.!  spectra  are  referred  to 
as  '>iiyn-rcL'o.,".ut!oji  i'pcjim.  for  spectra  computed  from  1 ,i'74-second  records, 
this  tiigli  resolutic.  i ,s  approximately  ().('(>1  hcrtc. 

The  program  used  to  com)nite  these  spectra  i .s  the  same  program  used  at 
Cf.llC  for  tile  analysi.s  of  field  data.  Many  spectra  of  field  data  computed 
with  this  program  and  suminariced  by  groujung  11  adjacent  sjicctral  iieriods 
are  given  h>-  Tliomiison  (1974), 

flic  effect  Oi  spillover  in  the  spectrum  is  shown  in  tiie  figures  of 
Appendix  D.  liach  spectrum  resulted  from  eomiiining  only  three  sinusoids; 
however,  energ)'  contributions  appear  at  from  fiv<'  to  nine  adjacent  spec- 
tral period'  . 


23 


The  poi'Lods  as.sij;ned  to  tlic  sinusoids  giving  rise  to  simulated  obsei’- 
vations  1,  2,  3,  and  6,  differed  by  exactly  tlirce  spectral  periods  (see 
Table  3).  The  spectra  for  tliese  four  observations  (Tigs.  I)-1,  b-d,  b-3, 
and  D-o)  e,\hii)it  three  maxima  separated  i'y  two  minima.  The  energy  at 
tliese  two  minima  may  be  interpreted  as  due  to  spillover  since  no  sinu- 
soids were  combined  with  the  corres]iondi ng  periods. 

The  high-resolution  s]iectra  from  field  wave  observations  are  dis- 
cussed later;  however,  the  use  of  the  minima  in  tliese  siiectra  to  esti- 
mate spillover  and  noise  is  discussed  liere.  Tlie  average  energy  at  tiie 
minima  between  25  and  7 seconds  in  the  liigli- resolution  spectra  of  the 
field  observations  was  used  as  a measure  of  spillover  and  noise.  Only 
spectral  periods  displaying  an  energy  content  at  least  twice  this  "back- 
ground" energy  were  interjireted  as  possibly  arising  from  physical  wave 
comiionents  in  the  wave  field. 

The  average  directions  resulting  from  the  10  tliree-gage  arrays  are 
given  in  i'alile  4 (last  column).  Only  results  of  comjnitat ions  at  the 
spectral  iieriods  closest  to  tlie  assigneil  periods  are  sliown.  The  table 
shows  that  for  wave  trains  1,  2,  and  5 the  computed  directions  for  tlie 
10  arrays  agree  with  the  input  directions  to  within  1°. 

The  directional  results  for  wave  train  4 are  correct  only  for  spec- 
tral period  62.  The  main  difference  between  this  train  and  wave  train  2 
is  the  narrower  spectral  width.  Wave  train  2 gave  the  correct  directions 
for  all  spectral  periods;  wave  train  4 did  not. 

Tor  simulated  wave  trains  5 to  8,  tlie  average  directions  seem  meiuiing 
less.  To  determine  whether  these  poor  results  were  due  to  n programing 
deficiency,  another  eight  sets  of  simulated  records  were  generated,  inter 
changing  periods  and  directions.  The  computer  output  for  the  simulated 
observations  is  in  Apiiendi.x  I).  This  appendix  and  Table  .5  should  be 
referenced  in  the  following  discussion  of  additional  simulated  wave 
trains . 


The  directional  results  for  the  sinusoids  with  periods  clastered 
around  1(>  seconds  were  of  the  same  tpiality  regardless  of  the  assigned 
direction.  However,  the  directional  results  for  the  S-second  sinusoids 
indicate  that  the  cajiability  to  sense  tlie  correct  direction  for  these 
shorter  waves  depends  on  the  orientation  of  the  three-gage  array  rela- 
tive to  the  direction  of  ]i rojiagat i on  of  tlie  incoming  wave.  The  result- 
ing directions,  which  differed  by  less  than  87°  from  the  assigned  direc- 
tions, are  given  in  Table  5.  The  top  of  each  column  in  tlie  table  shows 
tlie  shajic  and  orientation  of  tlic  array.  These  results  are  not  surprising 
since  tlie  effective  gage  separation  for  the  different  gage  pairs  varies 
with  orientation  relative  to  direction  of  wave  projiagat  i on . fable  5 
also  shows  that  the  more  nearly  equilateral  arrays  have  wider  direction 
il  i sccniab  i 1 i ty . The  design  considerations  for  the  array  indicated  an 
effectiveness  for  wave  periods  between  2.5  and  7 seconds  (Borgman  .and 
I’anickor,  197(1). 


l';iblc  4.  (A)iiiputational  results  at  closest  spectral  frei.|Ueticies  for  simulated  wave  trams. 


Wave 

Closest 

Spectral 

Amplitude 

Direction 

train 

spectral 

period 

frequency 
(1/1024  hz) 

U) 

(cm) 

L-  - 

1 

61 

16.79 

5.62 

I 20  S. 

64 

1 6.00 

7.55 

2(1  S. 

67 

15.28 

i 

5.70 

20  S, 

2 

61 

1 6.79 

5.62 

21  S. 

64 

16.00 

7.55 

15  S. 

67 

15.28 

5.70 

20  N. 

3 

62 

16.52 

5.52 

20  S. 

64 

16.00 

5.35 

20  S. 

1 

65 

15.75 

1.19 

20  S. 

4 

62 

16.52 

5.52 

21  S. 

64 

1 6.00 

5.35 

25  S. 

65 

15.75 

1.19 

36  N. 

5 

125 

8.19 

5.62 

142  N. 

128 

8.00 

7.55 

142  N. 

131 

7.82 

5.70 

j 

142  N. 

6 

125 

1 8.19 

5.62 

142  N. 

1 28 

8.00 

7.55 

78  N. 

1 

131 

7.82 

5.70 

162  S. 

7 

126 

8.13 

5.52 

142  N. 

1 28 

8.00 

5.35 

142  N. 

129 

7.94 

1.19 

142  N. 

8 

1 26 

1 

8.13 

5.52 

142  N. 

128 

8.00 

7.14 

76  N. 

1 29 

7.94 

5.83 

172  S. 

rill'  lonsisti'iit  .kkI  I'xait  ri'iinory  of'  assijjiiotl  directions  achieved 
for  the  simulated  1(>- record  wave  trains  is  in  part  due  to  the  use  of 
the  high  coni'ut  at  i ona  1 resolution  available  in  a large  computer.  Use 
of  less  exact  data  as  a\ailal>le  from  recoriling  instruments  is  ex])ected 
to  result  in  a less  consistent  and  accurate  recovery  of  the  true  direc- 
tion. Simulated  observations  1 to  4 were  rerun  truncating  the  computed 
profile  values  to  three  digits  as  is  commonly  available  from  recording 
systems.  The  effect  of  this  truncation  is  estimated  to  have  introduced 
an  error  of  the  order  of  Ml.ld?  centimeter  (0.05  inch)  in  the  instanta- 
neous values  of  the  profiles.  No  appreciable  differences  in  computed 
directions  resulted  l\v  this  truncation. 

3.  Identification  of  Wave  Trains  from  the  Hi  gli- keso  1 ut  i on  Spectrum. 

A wave  train  in  the  real  ocean  is  conceivably  made  up  of  several  wave 
components  with  nearby  periods  propagating  in  ap]iroxi mately  the  same  direc- 
tion. Simplistic  idealizations  of  such  wave  trains  are  exemplified  by 
simulated  observations  1,  3,  5,  and  7.  For  wave  trains  1 and  3,  the 
long-crested  wave  model  gave  the  correct  direction  of  wave  propagation 
for  all  ID  combinations  (within  5°),  not  only  at  those  spectral  (leriods 
closest  to  tlie  periods  of  the  sinusoids  comlMiied,  but  at  several  adjacent 
ones  on  either  side  of  these  periods  (see  Aiip.  I),  Figs.  I'-l  and  D-3).  At 
some  of  these  adjacent  spectral  periods,  the  contiibution  to  the  energy 
was  several  times  tlie  background  level  and  nearly  the  same  at  the  five 
gage  locations.  Thus,  it  can  be  assumed  that  a wave  train  in  the  ocean 
will  give  rise  to  a number  of  adjacent  spectral  periods  in  the  high- 
resolution  energv'  spectrum  with  energy  content  several  times  the  back- 
ground level.  This  background  level  can  be  estimated  by  inspection  of 
the  minima  in  the  energy  spectrum,  as  discussed  previously.  A criterion 
for  what  energy  level  will  be  considered  "h  i gh -energ>’  content"  can  be 
set,  and  groups  of  adjacent  spectral  periods  in  the  spectrum  with  high- 
energy  content  identified,  ibese  groups  imiy  each  be  assumed  to  arise 
from  the  presence  of  a wave  train  in  the  field  with  a mean  wave  jieriod 
within  the  range  of  spectral  periods  iti  the  group  (a  wave  packet).  The 
number  of  adjacent  spectral  [>eriods  in  each  grouj'  will  lu'  used  as  a 
measure  of  the  width  of  the  energy  peak  in  the  spectrum  and  indicates 
the  spread  in  periods  of  the  wave  train.  Ihe  spread  in  computed  direc- 
tions at  adjacent  spectral  periods  in  a group  is  ,in  indication  of  the 
degree  of  directional  organization  in  the  wave  tr.iin.  I.irge  spread  in 
directional  results  may  indicate  the  possibility  th.it  crossing  w.ive 
trains  with  nearly  tlie  same  period  are  present.  As  results  for  simulated 
wave  train  4 indicate,  the  long-crested  model  based  on  the  assumption  of 
a single  wave  train  at  each  frequency  is  not  suitable  for  a determi  n.it  i on 
of  wave  direction  in  such  cases.  Multiple  wave  tr.iins  at  a single  fre- 
quency may  result  from  refraction  around  a shoal  or  island  or  from  reflec- 
tion by  a vertical  wall. 

4 . Spectra  and  Direction  of  Wave  Propagation  for  Field  Data. 

The  energy  and  direction  of  wave  propagation  at  each  spectral  period 
wore  computed  for  44  field  observations  where  the  average  uncompensated 
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significant  wave  height  was  over  2 feet  and  the  discrepancy  of  individual 
standard  deviations  from  their  mean  was  3 percent  or  less.  Plots  of  the 
high-resolution  spectra  for  these  observations  are  in  Appendix  E.  The 
vertical  lines  represent  the  energy  contribution  at  each  spectral  period. 
The  background  level  for  each  observation  was  estimated  from  the  minima 
between  25  and  7 seconds  in  the  spectra.  Spectral  periods  in  this  range 
with  energy  content  above  twice  the  estimated  background  energy  were 
identified.  Contributions  to  the  energy  satisfying  this  criterion  at 
adjacent  spectral  periods  were  considered  as  arising  from  the  same  wave 
train.  The  nunber  of  adjacent  periods  in  each  train  was  used  as  a measure 
of  the  speatval  width  of  the  train.  The  energy  had  to  be  above  the  chosen 
level  at  all  five  gages  for  the  spectral  period  to  be  included  in  the 
group.  The  spectral  period  among  these  showing  maximum  energy  was  taken 
as  the  "period"  of  the  wave  train. 

Directions  were  computed  at  all  the  spectral  periods  in  each  train 
for  the  10  arrays.  The  total  spread  among  these  directions  was  found,  and 
an  average  total  spread  was  computed  for  the  trains  having  the  same  spec- 
tral width.  The  same  was  done  for  the  computed  directional  spread  at  the 
period  of  the  train. 

Twenty- five  percent  of  the  identified  wave  trains  had  total  direc- 
tional (computed)  spreads  above  100°  and  were  not  considered  further. 

For  89  percent  of  the  discarded  trains,  the  period  of  the  train  was  under 
9.4  seconds.  Thus,  all  trains  with  periods  under  9.4  seconds  were  dis- 
carded. 

Results  for  the  different  spectral  widths  for  the  trains  retained 
(280)  are  shown  in  Table  6.  The  second  column  in  the  table  gives  the 
average  total  directional  spread  for  the  corresponding  spectral  width; 
the  third  column  gives  the  average  directional  spread  at  the  period  of 
the  wave  train.  The  last  column  gives  the  number  of  wave  trains  having 
the  spectral  width  in  the  first  column. 

These  results  indicate  that  the  total  directional  spread  increased 
with  frequency  width.  Narrow  peaks  consisting  of  from  one  to  three 
spectral  periods  are  most  frequent,  and  the  spread  in  the  direction  at 
the  period  of  the  train  remains  relatively  constant.  Since  the  average 
spread  for  narrow-banded  trains  (width  <0.003  hertz)  is  21.8°,  it  is 
expected  that  three-gage  arrays  cannot  yield  directional  results  to  any 
better  accuracy.  Tlie  mathematical  exercise  in  Appendix  C shows  that 
array  dimension  is  a limiting  factor  as  to  the  shortest  period  for  which 
some  directional  discrimination  may  be  expected.  An  important  factor  in 
the  validity  of  the  directional  result  is  the  spectral  structure  of  the 
wave  train  involved.  Only  in  very  special  circumstances  will  the  quanti- 
ties involved  in  equation  C-10  (App.  C)  combine  to  give  better  results. 

There  are  various  possible  explanations  for  the  large  spreads  observed 
in  the  directional  results  from  field  records.  For  the  long-crested  wave 
model  to  be  strictly  applicable,  it  is  important  that: 
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Tabic  6.  Average  spread  in  euinpiited  direetioiis  lur  280  vva\  tiams  icleiuitied 
in  the  high-resolution  spectra  o(  4*1  field  wave  obs  -i \ a tions. 


(a)  Tlie  phase  differences  be  known  accurately  or  that  the 
probable  error  in  their  computed  values  be  known. 

0^)  The  wave  crests  over  the  array  site  be  long  and 
straight;  thus,  the  waves  must  not  have  undergone  appreciable 
refraction. 


(c)  'flic  sea  surface  be  stationary  in  time  for  the  dura- 
tion of  the  record  and  in  space  over  the  span  of  tlie  array. 

Tfie  mathematical  exercise  in  Appendix  D indicates  that  the  analysis 
yields  accurate  phase  differences  only  for  strictly  monochromatic  condi- 
tions. When  this  is  not  the  case,  no  accurate  estimate  of  the  error 
involved  in  the  computation  of  direction  can  be  given.  This  inability 
is  inherent  to  the  computational  procedure  and  cannot  be  resolved. 

Waves  with  periods  over  8 seconds  have  been  and  are  undergoing 
refraction  at  the  site  of  the  array.  Therefore,  tlie  wave  crests  are 
not  exactly  straight.  For  the  longer  waves,  with  wavelengths  at  tlie 
array  site  several  times  the  gage  separations,  the  curvature  will  not 
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introJuce  much  error.  This  will  not  be  the  case  for  the  shorter  waves 
and  orientation  of  the  array  becomes  important. 

Because  of  refraction,  the  curvature  of  a wave  train  changes,  perhaps 
only  slightly,  while  propagating  over  the  array.  This  change  will  intro- 
duce differences  in  the  direction  at  each  gage  and  differences  in  the 
direction  sensed  by  different  gage  pairs,  causing  undetermined  additional 
errors  in  tlie  computation  of  direction.  To  determine  the  magnitude  of 
these  errors,  two  additional  sets  of  simulated  wave  records  were  gene- 
rated. The  periods  of  the  sinusoids  combined  were  those  for  simulated 
observation  3;  the  directions  assigned  were  spread  within  a 10°  arc  for 
the  first  set  and  20°  for  the  second.  The  last  two  computer  outputs  in 
Appendi.x  D show  that  spreads  of  the  order  of  lb°  and  32°,  respectively, 
resulted  in  computed  directions. 

A stationarv'  condition  in  time  is  usually  assumed  when  developing 
wave  directional  models.  Indications  are  that  this  is  not  strictly 
applicable  at  all  times. 

Tlic  three  factors  discussed  above  are  sufficient  to  account  for  the 
inaccuracies  encountered  in  the  computations. 

5.  Conclusions. 


The  results  of  directional  computations,  for  both  simulated  and  field 
wave  data  records,  indicate  three-gage  arrays  have  some  capabilities  to 
determine  wave  direction  under  certain  conditions.  These  capabilities 
depend  on : 

(a)  The  dimension  of  the  array  and  the  water  depth  at  the 
site  which  place  a lower  limit  on  the  wave  period  for  which 
possibly  accurate  directions  may  be  computed. 

(b)  TIic  orientation  of  the  array  for  the  shorter  periods. 

(c)  llic  nature  of  the  wave  field;  directional  results  for 
wave  trains  with  a narrow  frequency  distribution  or  where  the 
computed  directions  differ  little  at  the  adjacent  spectral 
periods  might  be  meaningful. 

For  wave  trains  with  narrow  frequency  and  directional  width,  and 
period  above  10  seconds,  the  three-gage  arrays  at  Pt.  Mugu  yield  direc- 
tions to  ;in  estimated  accuracy  of  20°. 

At  the  Pt.  Mugu  site,  16-second  waves  may  approach  the  coastline  at 
angles  of  22°  or  less  from  the  normal.  Tlie  directional  information  pro- 
vided by  the  array  adds  little  to  this  and  seems  hardly  cost  effective. 
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APPENDIX  A 


DERIVATION  OF  THE  EXPRESSION  FOR  WAVE  DIRECTION 


Let  the  coordinates  of  nearby  gage  sites  be  i = 1,  N; 

with  N = number  of  gages.  The  water  surface  displacement  at  each  site 
due  to  the  passage  of  a sinusoidal  wave  of  frequency  o,  and  amplitude 
A,  traveling  in  direction  a (a^  in  Fig.  8),  is  given  by: 


A cos{k[x^  cos  a + y^  sin  a]  - 2not  - iji}  , 


(A-1) 


where  k = 2'tr/L  is  the  wave  number,  L the  wavelength,  and  (ji  the 
initial  phase  at  the  origin. 


Tlie  phase  difference,  • 

sinusoid  considered  is: 

between 

locations 

i 

and  j 

for  the 

Tims,  for  three  noncolinear  locations. 

(Xi  - Xj) 

s in 

a]  . 

(A-2) 

$12  = k[(xi  - X2) 

COS  ot  + 

(Xl  - X2) 

sin 

a]  > 

$13  = k[(xi  - X3) 

cos  Qt  + 

(Xl  - X3) 

sin 

a]  , 

(A- 3) 

suffice  for  a unique  solution  of 
sin  a terms  and  then  the  cos  a 

the  direction  a. 
terms,  to  obtain: 

El 

iminating 

first  the 

sin  a = 


cos  a = 


Cxi  - X3)  <l»i2  - (xi  - X2)  <}>13 

J<[(xi  - X3)(yj  - y^)  - (xj  - X2)(yi  - yj)] 

(Xl  - Yi)  <{’13  - (Xl  - X3)  *^12 

kf(xi  - X3)(yi  - yz)  - (x:  - X2)(Xl  - X3)] 


(A-4) 


Since  k is  always  positive,  consideration  need  only  be  given  to  the 
other  terms.  Letting  D stand  for  the  quantity  in  square  brackets,  the 
direction,  a,  for  D ^ 0 is  given  by: 


a = tan 


- 1 


[Cxi  - X3)  - (X;  - X2)  »I3]/D 

Kxi  - X2)  '5’13  - (Xl  - X3)  ^12]/^ 


(A-5) 


A unique  value  for  direction  can  be  obtained  by  considering  the  signs  of 
both  numerator  and  denominator.  The  quantity  D differs  from  zero  for 
all  nonlinear  arrays  as  shown  below. 

Let  x^  = yj  = 0,  D will  equal  zero  for  yz/^2  ~ X3/X3;  thus, 

X2 , y2  and  X3,  y3  will  be  on  a straight  line  with  slope  given  by 
this  ratio. 
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APPENDIX  B 


COMPUTER  OUTPUT  FOR  CROSS-SPECTRA  COMPUTATIONS 


Table.  Guide  to  computer  output. 


Line  ' Column  Explanation 

(from  top)  (from  left) 

1 Title,  plus  date  and  time  of  the  observation 

(day,  month,  year,  hour,  and  minutes) . 

2 to  4 Headings  for  columns.  The  numbers  separated  by 

dashes  in  the  fourth  line  give  the  numbers  of  the 
gages  in  the  array  (see  Fig.  3). 


Sequential  nunfcer  of  bands  0.0107  hertz  wide. 

Period  at  center  of  band  (seconds) . 

Percent  of  energy  in  each  band  for  gages  1 to  5. 
Normalized  to  the  energy  content  from  approxi- 
mately 30  to  3 seconds. 

Resulting  "representative"  direction  of  wave 
propagation  for  the  10  arrays  for  the  corre- 
sponding band  (degrees) . 
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APPENDIX  C 


FOURIER  COEFFICIENTS  FOR  A MIXTURE  OF  THREE  SINUSOIDS 


The  Fourier  Transform  of  the  function: 

f(t)  = A cos  (at  - (}i) 


(C-1) 


with 


2ir(m  + 6)  i t 1 
NAt ' 2 ’ 


(C-2) 


where  1 <<  m <<  N,  confuted  from  N values  of  f(t)  evaluated  at 
equal  increments  of  t;  At  is  given  by  the  set  of  coefficients  (Harris, 
1974) : 


- 


= 


2A  sin  7r6  cos  (ttA  - 4>) 
N 

2A  sin  7t<5  sin(7i6  - cb) 


N 


r 1 ^ 1 ] 

|tan[n(m  - m + 'S)/NJ  tan[n(m  + m + 6)/N]J 

r .1, 

[tan[n(m  - m + i5)/N]  tan[iT(m  + m + <5) /N]J 


m = 1.  2.  . . . 


(C-3) 


Harris  shows  that  for  values  of  m near  ra,  and  for  m far  removed 
from  1 and  N/2, 


an  = 


A sin  tt6  cos((j)  - 1:6) 

IT (m  - m + 6)  ’ 


bn  = 


. A sin  tt6  sin(<{i  - Tr5)  , ^ m 

. ra  = 1,  2,  . . . _ 

TT(m  - m + iS)  2 


(C-4) 


are  good  approximations  to  the  coefficients  in  equation  (C-3) . Equation 
(C-4)  shows  that  convergence  is  slow. 

For  the  special  case  6=0,  substitution  into  equations  (C-3)  or 
(C-4)  gives: 


and 


- 0 > m ^ in  , 


an  = bn  = indeterminate  for  m = in 
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Use  of  L'  Hospital  Rule  in  equation  (C-3)  when  6 ->■  0,  shows  that: 

a^  = A cos  (j)  , = A sin  ij>  , m = in  . (C-5) 

Application  of  the  cosine  bell  data  window, 

f(nAt)  = 7 j[l  - cos  ^]j  f(nAt)  , n = 1,  . . . N . (C-6) 

is  equivalent  to  replacing  the  original  sinusoid  f(nAt)  with  the  sum 
of  three  sinusoids  (Harris,  1974),  where 

re  A L)  f2Tr(m  + 5)n  ,1 

f(nAt)  = 2 2 cos  I-  - - <J> 

r2Tr(in  - 1 + 5)n  1 r2Tr(m  +1  + 6)  l| 

- L — s ^ ■ !J  ■ ^ ■ tJj  ■ 

For  6=0,  and  in  view  of  equation  (C-5),  the  Fourier  Transform  of  this 
modified  function  will  be  given  by: 

A A 

^m-1  ~ ^m+1  ~ ~ J ” ■ ^4  'J' 

and 

A A 

a^  = j cos  <ti  , b^  = - sin  4> 

with  a^  = = 0 for  all  other  values  of  m.  Thus,  energy  appears  at 

three  adjacent  m values. 


Harris  (1974)  shows  that  after  application  of  the  cosine  bell  data 
window,  the  approximate  values  of  the  coefficients  are  given  by: 


A sin  ti6  cos(<1)  - tt5) 

2Ti(m  - m + i5)  [(m  - m + 6)^  - 1] 


A sin  it6  sin(<Ji  - it6) 

2Tr(m  - m + 6)  [(m  - m + 6)^  - 1] 


m = 1 , 


(C-8) 


Thus,  convergence  increases  rapidly,  and  values  of  the  coefficients  for 
(m  - m)  ^ 3 may  be  disregarded. 

Equations  (C-3),  (C-4) , and  (C-8)  imply  that 


tan((()  - 7i6)  = ^ 
bffi 


(C-9) 
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Since  6 may  be  as  large  as  1/2,  phase  values  computed  from  the 
Finite  Fourier  Transform  may  be  in  error  by  as  much  as  90°. 

For  a simplistic  simulation  of  a wave  train,  it  is  sufficient  to  com- 
bine three  sinusoids  with  nearby  periods  propagating  in  the  same  direc- 
tion. Letting  equal  the  amplitudes  and  k^,  i - 1,2,3,  the  wave 

numbers,  the  Fourier  Transform  of  the  combination  is  given  by; 

^ ^ H sin  tt6-^  cos(<i>/  - ■nS>j) 

2Tr(m;  - m + - m + 6^)2  . 


ki  sin  sin(4>£  - 
i=l  - m + <S^)[{%  - m + 6^)2  _ 


Nearby  periods  are  attained  by  setting 


= (m^  - m)  < 3 


CC-10} 


(C-11) 


Let  the  coordinates  of  three  nearby  locations  be  (xj,yj) 5 j = 1,2,3. 
The  only  difference  among  the  Fourier  Transforms  (eq.  C-10)  arising  from 
wave  records  at  each  location  is  the  values  of  the  At  each  loca- 
tion j,  the  values  are: 


~ a + yj  sin  o) 


i = 1,2,3  , (C-12) 


where 


ij  = Xj  cos  a + yj  sin  a , j = 1,2,3  . 


(C-13) 


Since  the  three  sinusoids  are  assumed  to  have  nearby  periods,  let 


ki  = k2  = k3  = k . 


: !ius , 


Then,  for  the  wave  record  at  location  j : 


3 sin  it6^  cos(kn^-  - 
a-mj  = S *■ r 

i=l  2ir  + 6/3  [(A/  + - 1] 


3 sin  sin(knj  - 
” i = l 27f(Ai  + 6^3  [(A.^  + - 1] 


(c-143 
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l.ct ; 


and  sin  k J , is: 


+ sin  k 


and  s imi  lar  ly  for  b,.  r : 


b ■ = sin  k : 

■a 


- cos  k. 


• + 

- 1]  . 

i 

= 

, a f t e r 

col lec 

ting 

te 

(At 

s i n 

. . 2 

cos (it 

+ -I  • 

l) 

V ' 

1 

[1] 

A2 

s i n 

•I  ' 

cos  ( i 

+ * 

s) 

A 3 

s i n 

r, 

o — 

+ -6 

3) 

I3i 

(Ai 

s i n 

s i n ( : J 

+ r '■ 

1) 

1 

m 

A 'I 

s i n 

s i n ( ; , 

+ -f. 

.’1 

U1 

At 

s i n 

t 

s i n ( : 3 

+ -6 

tl 

13] 

|At 

s i n 

" ■ 1 

■/. 

0 

+ -6 

1^ 

t I 

m 

A. 

s i n 

cos  ( t 

+ •:  • 

l-'l 

At 

s i n 

COS  ( i 3 

+ - • 

13) 

I 

At 

s i n 

■’■■1 

s i n ( 1 1 

+ -•  ' 

c 1 

m 

A;' 

s i n 

s i n ( i , 

+ - • 

,0 

l-.'I 

At 

s i n 

" ' t 

s i n ( J 3 

+ 1 ■ 

t' 

for  every  narrow-banded  wave  train,  I | ■'  3 in  ei|uation  (Ci-lS), 
unprcdictal)  1 e terms  arc  introduced  in  the  ex[iressions  for  the  coefficients 


5I 


k 


wilich  make  the  ratio  h„/a,„  a poor  estimator  ot'  the  value  of  the  phase 
of  the  sinusoid. 


'llie  numerators  of  the  terms  inside  tlie  braces  in  tlic  alio\X'  cc|uations 
for  a,,-  and  b,„ : are  at  most  of  order  A^,- , i = 1,2,  a. 

Randomness  i:i  the  values  of  i = 1,2,3  and  in  tlic  phase  rela- 

tionshi[is  of  tile  three  sinusoids  at  tlic  origin  of  coordinates  might  pro- 
duce partial  cancellations  among  tlie  terms  inside  the  braces  to  reduce 
tlic  resulting  error. 

-\ssune  for  example:  = 0,  i = 1,2,5.  The  coefficients  reduce  to; 


.'1  • = 

~ cos  kfi  : 

sin  2t1  j A; 

+ 

sin  2-  ' . 

‘'^3 

sin  2 ■ 

3I 

1 

m 

12] 

|3] 

1 

+ sin  i;.-.- 

J 

h 

s i n - S 1 A 

. sin^-  V. 

■I- 3 

sin''  ‘ 

3} 

\ 

\ 

11) 

[2] 

13] 

1 * 

s i mi  1 arl>’ , for  b 

b,„  • = 

■ 

1 - , 
“Sin  k : 
J 

(-M 

sin  1 A; 

1 sin  2-‘- 

■'3 

sin  2tt6 

.)  1 

1 

[1] 

[2] 

13] 

- cos  k;j  ; 

I-'^l 

siir  A. 

, sin*^  - 

A3 

sin^  tt5 

3 1 

1 

(C-16) 

in 

[2) 

13] 

1 

f • 

Letting 

I-l  = 

Ai  sin  2-’ 5 

1 

) 

1,-2  = Aj  sin"  - 1 , 

Ml  = 

A.  sin  2"'. 

) 

M;.  = A;  sin-  ■■■■  , 

\l  = 

A 3 sin  2r  ■■ 

3 

J 

.N-  = A 3 sin-  ••■■.3  : 

(C  17) 

tlien : 

-dr-'  = 

1 cos  k-: 

2 

Fli 

. Ml  , Ni  ■ 
[2]  [3]_ 

+ sin  k _ r 

F-/ 

jl] 

''c 

" |1T 

+ 

13]_ 

f 

and 

i s i n ki. ; 

Ml  Ni  1 

. + -f  - 

[-’)  [3]_ 

- cos  k : • 

1,,. 

M-, 

' V] 

4- 

.N-, 

[3]_ 

(C-18) 

Assume 

furtlicr : 

Ai  ~ 

A^.  ~ A3  and 

1 m ; - m 1 = 

0 

for  i 

= 

2 and 

equa  1 

for  i 

= 1 and  5. 

Since  I'j  •'  1, 

tlie  terms 

m 

and 

[5]  are 

of 
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approximately  the  same  magnitude  but  of  opposite  signs.  Thus,  terms 
involving  the  products  [1]  [2]  and  [2]  F3]  tend  to  cancel.  Letting 

6 j = 6 o = 'S  2 = 5 ' 

h„-,-  sin  kfi  ; sin  2tjc  - 2 cos  k2  • sin^  ’t5 

• t'  c'  ^ 

~ cos  kP.,-  sin  2-mS  + 2 sin  kP.;  sin2  tt6  ’ 

Using  the  trigonometric  identities  for  the  double  arc,  this  expression 
reduces  to: 


— — = tan  (kP,-  - TTC)  . 


(C-19J 


Phase  differences  between  locations  will  be  approximately  correct. 
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Al’PHNDIX  1) 


SPECTRA  PLOTS  AND  COMPUTER  OUTPUT  FOR  SIMULATED  OBSERVATIONS 


F'igures  D- 1 to  D-8  sliow  high-resolution  spectra 
pressure  gages  1 to  5 at  Pt.  Mugu,  California. 


54 


r tr  W V • • ) 


/ 


Figure  D-b.  j 


60 


r 


ijHut*Tto  SE»  «>e«iansfsec)  oiocc^ioxseotcBEts)  iBBunuDiscc") 

16. M 21.  II. uJ 

11.  11.26 
11.26  .21.  n.us 


I -.-4 

i 


»l^»»9<^9»9<^9c^«)«c^r»«lwAl•-6Moatf>9Hlo••a>nJtf'o«%•«M«•eooooooooooc»o 

*r»» 

<^^f^9<^99»l^^»99«^«3<^^••M09(A9'^0•«■>(^lflOI^«fli>•*•0000000000000 

«»»»• 


•'90'<^9^<r<^<r9<^c««^»>Afin«*^•-•oo9tf'9Klo•>•^'Mtf^^«•ooooooeooooooooo 

»a^o<B99<r9'9»<B««4i99V>Ai*^«««>*o4ir9K6o.«r'Ot^*u^oooooooooooooeoo 


»e*«o<wo^** 


SI**UU*TID  XIUD  set  >>FB10DS(8EC)  Oinec^ IOxsFOECRECS)  o-'LltuOtlCCX) 

I*!***  21a  Hall 
l*a«T  21a  lta2« 
>Saf2  21a  lla*> 


ti  v / 


I 


• r.  ^ « iA  4^  <W  <V 


ru 


Q _A 
QC  «>i  9 

•*  ac  • 

K ,M  M 


« ««l/>lAtrir<riA(Atr4A  9 «miAiAiAAlAlAlA*^A«AifWIWM<W«WI^AllWMI«W(WM<MAJMM<WM«B»««*««««>« 


«0  3 
O 

UlA  9»0»»9Atfa2»«Ar>'*^lA4^0*0»*M«*00»«0»Afy*lAlA<VIV*«*«««00000000000 

orat  ■•»•»•••••••••••••• 


« A <A  A « A 


9 « A « 


9Al(A0fA9^««0«»00>»«aA»te 

~ ^ • 

B • • 


a 


9 «4B»^B^4B9lA«A4«««««AO»«.|A«*«««^MO««OOMO»B>-lAl4««*0»*B>A|B»«a»A«M^A|M>««««»-* 

f\l  «-•  M •.  M M M M •.  .M  M ••  M 


Ai  tBBBtBBBBBB«*^«>**«»«>«*> 

B BBBBBBB 


ooooooooooooooooooo**«a*b^«999oooooc»oocboooooooo»ooc»o 

00000000000000000000#<B^490000000000e000000000004BO 


lA  o m « ••  o 

fA  9 «• 


OOOOOOOOOOOOC-OOOOOOOB 


10000000004 


« O «A  V «•  O 


>000004 


>00000000000 

>00000000000 


oeooooooooooeoooooo«"*MOB^*M(WM900oooooooooooooooooooo 

0O0000000eoO00O00OOO99lAB»  9000000000000000000000000 


ir  o ^o  ^ o 


O M 
MB  •» 

& W 


09«^«nAIA4*«**e0O 


65 


SI«iL»TEO  “I«En  SE*  PF»lonStsEn  CI«rCElP*.S(DEC«E£S)  »“Pl1TI.B£S[C») 

I*."*  ?1.  II. «s 

1S.«T  li.  15. 


- fV  rw  ' 
• «w  f\<  <V  •*  t 

• • • 


>00000000 


><rC^O99»«t40«O»'Ai«/VOO«>•»>■lV«>A«AJ9O«^ 

' — — ••• 


l^t^f^MMOOOOOOOOOOOOOO 


rw'\jr\j<\jr\j'\ir\jrvjrxj'\i'\jK5ti^  ^ :>r^oo^»^09rv«~*« 


^ ■O  o O a ^ t r\4  ^ 


*0000000 • 


a»otfaaacr»9o< 


’ Cl  «/  A a o 


4 9o»»»<r»<r»<r»9<c>e«i^«>a»>MO»* 

t ► r-  ►•  r- 


>0000000000000  c> 


a — • O O «»  <c  9 >^fS<«K\C>0«»< 
- — - — — — • I 


— — r,,^  11—  111-11  - .-  i.i  i.  .1 


a a a a a a a a amtr«/>ir>tf*i/>v>t/««r‘i/t«««  « 


66 


• i"ut»'io  «i»io  M*  Bi»ecTiom(i)i*«K»)  *«»ufuot«c») 

• .It  •••.  II-*! 

t.M  •40.  It. 14 


67 


■ io"Sf  lec)  of»ecTio*«s(ot  «•€»•)  ^w^lituocscc*) 


72 


9 


>*'»'9aa99'iaaa9-«' 


oeoo»  oi^«ma«^«w*w<ma««w»>mooooo  oooooooo* 


A^r>~ckMO-*>^o^r»  r^  9 a 9 o 


^ »«.«  »>r>  9 w ^ (»^*MAjMi«>-*eo^ 


■ . ^fV*VA*A>A#<V*V^rwA<«A«Ap^0«>AlA»A^V 


T10N9(0tei>ECS)  *>^LITuOe»(C>i) 


r 


r 


« o o o o 


M ^ ^ ^ • 


M «W 
O •-« 


X »-  • 


^1 

:s 


• rV«/'«34>9t 
»•>  4)  ^ 4 r<-  y>  ^ . 

» • • 


‘•fKjH%av'(04>4 


• ^ ^ . 


xT  •“  — c*  4 c 


M’PliNDIX  i; 


lilGll-RliSOUrnON  SPIICTUA  FOU  Fll-I.l)  WAVF  [)A1A 


Figures  F-1  to  li-44  show  liigh-rcsolution  spectra  for 
pressure  gages  1 to  5 at  Pt . Mugu,  California.  Date 
ami  significant  wave  height  are  indicated  for  each 
figure. 


Energy  Density  (cm  / Hz) 


Period  ( s ) 
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0.06  0.10  O.H  0.18 

F requency  ( Hz ) 


ijMirc 


S i j;n  i fi  c.'int  w.-ive  heiszlu  .‘U.O  centimercrs  ( 
.’It  April  I'.f’i),  1/JS  lioui'S. 


0.06  O.IO  O.U  O.IS 

Frequency  ( Hz ) 


Figure  Ii-6.  Signi  f ic;uit  wave  height  100.7  centimeters  (3.4  Feet), 
21  April  1970,  0929  hours. 
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liguro  H-l).  Significant  wave  lieight  85,1  centimeters  (2.8  feet), 
21  .April  11)70,  2129  hours. 
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Period  ( s ) 


Frequency  ( ^ z ) 


ii'.uro  li-lii. 


S i j;n  i fi  cant  wave  hr'ii’Ju  SS.(>  cent  i mi-rt- 1- 
11  Illlie  r.)''n,  hours. 
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rijMiro  1-20.  Significant  wave  hcigtu  72. S centimeters  (2.1  feet), 
24  June  1970,  1807  lioiirs. 
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Fi}',uro  l'-23.  S i j’.n  i f i c;mt  wave  hciqht  6(i.7  centimeters  (2.2  feet), 
23  June  1070,  1808  hours. 
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Fijjure  [;-41.  Significant  wave  heiglit  91.2  centimeters  (3.0  feet), 

1644  hours. 


0.06  O.tO  0.14  0.18 
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Figure  E-44.  Significant  wave  height  G6.9  centimeters  (2.2  feet), 
18  Occember  1970,  0745  hours. 
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